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I derive uni-directional wave equations for fields propagating in materials with both electric and 
magnetic dispersion and nonlinearity. The derivation imposes no conditions on the pulse profile 
except that the material modulates the propagation only slowly: i.e. that loss, dispersion, and 
nonlinearity have only a small effect over the scale of a wavelength. It also allows a direct term-to- 
term comparison of the exact bi-directional theory with its approximate uni-directional counterpart. 
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I. INTRODUCTION 

In the past few years, composite materials ("metama- 
terials" ) that demonstrate both an electric and magnetic 
response have been the subject of both experimental and 
theoretical investigation. Often the motivation for this 
research is the potential for exotic applications: for ex- 
ample, superresolution [l| , or the possibility of "trapped 
rainbow" light storage Despite these interesting pos- 
sibilities, there is also the more basic need for efficient 
methods for propagating optical pulses in such metama- 
terials, in particular in one-dimensional (ID) waveguide 
geometries. Indeed, methods for doing so have already 
led to interesting predictions (see, e.g.Js'--^). However, 
these methods, and earlier ones, (e.g. [6|-t9i]) tend to rely 
on mechanisms such as the introduction of a co-moving 
frame, and assumptions that the pulse profile has neg- 
ligible second-order temporal or spatial derivatives. As- 
suming second-order derivatives are small may well be 
reasonable, but it means that the pulse profile must re- 
main well behaved. This approximation therefore might 
well become poorly controlled [13] , particularly for ul- 
trashort or otherwise ultrawideband pulses, or exotic or 
extreme material parameters. Ideally we would prefer 
to make approximations based solely on the material pa- 
rameters of our device, so as to avoid making assumptions 
about the state of an ever-changing propagating pulse. 

Here I derive ID wave equations for a waveguide with 
both electric and magnetic dispersion, and electric and 
magnetic nonlinearity. I use the directional fields ap- 
proach [ll|, , which allows us to directly write down 
a first-order wave equation for pulse propagation with- 
out complicated derivation or approximation. We simply 
look at the coupled forward and backward wave equations 
that are a direct re-expression of Maxwell's curl equa- 
tions, and substitute in the appropriate dispersion and 
nonlinearity. I also show separate examples for second- 
and third-order nonlinearities in both electric and mag- 
netic responses, although the effects can be combined if 
desired. Note that these directional fields are applica- 
ble to more than just pulse propagation, as they have 



been used to simplify Poynting-vector-based approaches 
to electromagnetic continuity equations [l3l |. 

The derivation makes only a single, well-defined 
approximation to reduce the bi-directional forward- 
backward coupled model down to a single first-order wave 
equation - that of assuming small changes over the scale 
of a wavelength. This approximation is remarkably ro- 
bust for all physically realistic parameter values - see 
for an analysis focused on nonlinear effects; more general 
considerations have been dealt with in terms of factorized 
wave equations jl5| . The resulting wave equation retains 
all the usual intuitive and analytical simplicity of ordi- 
nary wave propagation equations, unlike the computa- 
tionally intensive approach of a direct numerical solution 
of Maxwell's equations (see, e.g., [tI. [l6l - [20l |V 

This paper is structured as follows: Directional fields 
and their re-expression of the Maxwell curl equations is 
outlined in Section [III followed by the reduction of the 
bi-directional wave equation into a uni-directional form 
in Section IIIII Section IIVI shows wave equations for a 
doubly-nonlinear third-order nonlinearity material, and 
Section |V] does the same for a second-order case. In Sec- 
tion |VT] propagation under the influence of typical meta- 
material responses is discussed, and I conclude in Section 

rvm 



II. DIRECTIONAL FIELDS 

The directional fields approach ill] allows us to write 
down wave equations for hybrid electromagnetic fields 
G^. Note that here I define with an alternate (and 
more sensible) sign convention than previously. Further, 
I also allow for more general types of polarization and 
magnetization in such a way as to provide a simpler pre- 
sentation. For propagation along the unit vector u, the 
propagation (curl) equation for G* is written in the fre- 
quency domain as^ 

V X G^ = ±iWar/3rU X G^ ± lUJ^rVi X Pc -I- lUJar^^oM.c 

(1) 



'Electronic address: " Dr . Paul. Kinsler ©physics . org| 



^ See derivation in Appendix 1X1 
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with 



G=^(a;) = a^(w)E(w) ± /3^(w)H(w) 



X u. 



(2) 



Here the electric response of the material is encoded in 
two parts: a spatially invariant linear response compo- 
nent ar, and the remaining contributions (of any type) in 
Pc. Similarly, the magnetic response is divided up in the 
same way between f5r and /ioMc. Generally we will put 
the entire non-lossy linear response of the medium (i.e. 
the dispersion) into the reference parameters and f3r, 
although it may also be convenient to specify only that 
the product a^/Jr is real (cf. [2l|). All the nonlinear re- 
sponses and other complications ( "corrections" ) , such as 
spatial variations in the material parameters, remain in 
Pc and Mc. As an example, in this approach was 
applied to second-harmonic generation in a periodically 
poled dielectric crystal. The time derivatives of these 
corrections Pc and Mc correspond to bound electric and 
magnetic currents respectively [Ts^. These Pc and Mc 
are functions of both fields E and H, i.e. Pc = Pc(E, H) 
and Mc = Mc(E,H). If we choose instead to have 
frequency- independent and then the remaining 
linear response can simply be included in Pc and Mc; 
in this case the "weak loss and nonlinearity" condition I 
use later to decouple forward and backward fields would 
then need to be broadened to include weak dispersion as 
well (also see [ll| for more discussion). However, neither 
version imposes any requirements on the pulse profile. 

It is useful to give a simple example of the directional 
fields to provide some insight into their nature. In the 
pure transverse plane-polarized case, with fields propa- 
gating along the z direction, and frequency- independent 
(material parameters) permitivitty Cr and permeability 
Hr, we can write 



(3) 
(4) 



where this simple definition matches the original pro- 
posal of Fleck [2I] . 

It is worth considering how refiections arise in this 
picture based on spatially invariant reference parame- 
ters augmented by corrections terms. Leaving aside for 
now the distinctions between spatially propagated fields 
and temporally propagated ones (see the discussion in 
[isj). transition to a new media can be handled in two 
ways. First, we could map the existing fields onto 
new ones G^ based on new reference parameters am 
and (3rn- Here a pure G'*' field would separate into two 
pieces, one a forward propagating G^, and the other a 
"reflected" backward propagating G~ . Second, we might 
retain the existing reference parameters, and have mod- 
ified corrections Pc and M^. These altered correction 
terms then couple the forward and backward directed 
fields, inducing the necessary reflection in G~; although 
as a side-effect of our now no longer optimal and /3r, 
the forward evolving field is made up of coupled G^ and 
G~ components 



A. Material response 

We define the electric and magnetic material response 
in the frequency domain, as it greatly simplifies the de- 
scription of the linear components. Let us chose a ref- 
erence behaviour given by 6^(0;), /^^(a;), and use them 
to define reference parameters 0:^(0;) = ^/e^^uj) and 
/3r(w) = fJ,r{<^)- Notc that these are allowed to have 
a frequency dependence [llj], and that arPr is just the 
reciprocal of the (reference) speed of light in the medium 
(i.e. Tir/c). We therefore have that the electric displace- 
ment and magnetic fields are 

D(tj) = eoE(tj) 4- P(w) = e^(a;)E(a;) +Pc(a;), 

(5) 

B{w) = /xoH(w) + M(w) = /Ltr(w)H(tj) + ^oMc(a;). 

(6) 

To give a specific example, we can define frequency- 
dependent loss and dispersive corrections by Ke{uj) and 
K^(a;), along with (e.g.) independent third-order nonlin- 
earities Xe,Xfj, to both the material responses; although 
any appropriate expression can be used - even magneto- 
electric or other types. Thus we can write the frequency 
domain expressions 

Pc{uj)=alK,-E + eoJ[XeE^{t)]*-E (7) 
MoMc(a;) = p^.K^U + /xq J [xt.H^{t)] * H, (8) 

where ?[...] takes the Fourier transform (which is nec- 
essary because nonlinear effects are defined in the time 
domain as powers of the field) and ★ denotes a convolu- 
tion [i.e., a-kb = J a{uj)b{u! — uj')du!']. If the nonlinearity 
is time dependent, then the simple XeE^ type terms can 
be replaced with the appropriate convolution. In general, 
it is best to pick ar,(3r subject to the condition that the 
sizes of Pc and Mc are minimised. 

In a double-negative material (with both €,fi<0) we 
would get imaginary ar,/3r, changing the complex phase 
of G^(a;) away from that given by the original E and H. 
Since this is in the frequency domain, it converts into a 
phase shift in the time domain, so although imaginary- 
valued ar,f3r might seem inconvenient, it does not give 
unphysical results. 



III. WAVE EQUATIONS 

Starting with the vectorial curl equation ([T]), I first take 
the ID, but bi-directional, limit, and describe the approx- 
imation necessary to produce a simpler uni-directional 
form. After this, I discuss how the common transforma- 
tions used in optical wave equations can be applied in 
this context. All equations and field quantities are in the 
frequency (a;) domain, unless explicitly noted otherwise. 
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A. Bi-directional wave equations 

Here we set u along the z axis without loss of gener- 
ality, and consider just an x polarized wave (i.e., con- 
sisting of Ex,Hy). This means we use the y component 
of eqn. ^ with dz = d/dz, so that the wave equations 
for the full spectrum fields G^(a;), coupled by corrections 



Pc.x{^) and My 



Mc.y{uj) are 



iLoarMy. (9) 



Following the detailed discussion in [lll^, we say that 
this wave equation propagates ( "steps" ) the fields forward 
along the z direction using oppositely directed fields G+ 
and G~ . These fields can be written as functions of ei- 
ther time or frequency, and pulses they describe therefore 
evolve (travel) forward or backward in time. 

Consider the example case with parameters n^iuj) and 
K^(uj), and Xe, Xt^ in eqns. ([7]) and (|H]). Defining kr{u}) = 
U!ar{u!)l3r{uj), wc get 

dzG"^ = ±ikrG^ ± — ^ — (G+ + G^ ) H ^ (G+ - ^ 



±^Jj[xe£;^]*(G^ + G-) 
-l-^^j[x,i?^]^(Gj-G, 



(10) 



Note that even for a frequency-independent choice of the 
reference parameters Ur and /3r, the reference wave vec- 
tor kr retains a (linear) frequency dependence. Also, the 
dispersion and/or loss parameters directly re- 

lated to e and /z respectively, and not to a refractive index 
n or wavevector k. This is why there is a factor of 1/2 
associated with their appearance in eqn. PH)) and subse- 
quently. 

If written in the time domain, these wave equations 
are seen to propagate the full temporal history of a field 
forward in space. There, the reference propagation given 
by ±ikrG^ becomes a convolution if kr retains a nontriv- 
ial frequency dependence. However, if we expand kr{Lu) 
around a central frequency uji in powers of cj — cji, we can 
instead convert it (in the time domain) into a Taylor se- 
ries in time derivatives, which is a popular alternative to 
the frequency domain form used here. However, if imple- 
menting a split-step Fourier method of solving these wave 
equations, dispersion is applied in the frequency domain, 
so that in general such an expansion is an unnecessary 
complication. 



of one wavelength - or, if you prefer, over time inter- 
vals of one optical period. This translates into a weak 
loss and nonlinearity assumption; and if the correction 
terms Pc and Mc include dispersion, a weak dispersion 
assumption is also made. These are rarely very stringent 
approximations. If |Pc| << |D| and |/^oMc| << |B |, 
then a forward has minimal co-propagating G~ pj| . 
Further, the forward field has a wave vector kr evolving 
as exp(-l-zfcrz), but any generated backward component 
will evolve as exp(— ifc^z). This gives a very rapid rela- 
tive oscillation exp{—2ikrz), which will quickly average 
to zero. Nevertheless, although achievable optical non- 
linearity coefficients fall well within this approximation, 
care may need to be taken with the dispersion, particu- 
larly if near a band edge or in the vicinity of a narrow 
resonance. 

A directly comparable approximation is treated ex- 
haustively in [15|, where although applied to bi- 
directional factorizations of the second-order wave equa- 
tion, the physical considerations are exactly the same: 
Deviations from the reference behaviour over a propaga- 
tion distance of one wavelength should be small. Note 
that the slow evolution approximation applied here is 
not the same as other "slowly varying" types of approxi- 
mation [e.g., the slowly varying envelope approximation 
(SVEA)] - although the physical motivation is similar, 
the approach used here is far less restrictive. 

After we apply this weak correction or "slow evolution" 
approximation, we set the initial value of G~ = 0, and 
can be sure that it will stay negligible. Thus eqn. (|9]) for 
the full spectrum, forward directed field G^{Ld) can be 
written as 

dzG^ ~ +tUJar(3rGl^ +tUj/3rPx + tOJarfioMy. (11) 

Alternatively, we can scale the G^ field so that it has 
the same units and scaling as the electric field, using 
F+{uj) = G+{uj)/2ariuj). This gives 

dzF+ +iwarPrF+ +^'^P^ +i'^ti^M+. (12) 

Note that E+ ^ G+ I2ar = F+ and H+ = G+/2/3,r = 
F+a,//?r, since G~ =0. In either version of these 
uni-directional equations, P+ = P{E^ ,H^) and M+ = 
M{H^ , E^) - the uni-directional (residual) polarization 
P^ and (residual) magnetization AI^ should not be writ- 
ten as functions of the total fields Ex and Hy . 



B. Uni-directional approximation 

Now we apply the approximation: that the effect of 
any correction terms is small over propagation distances 



Section III 



C. Modifications 

Either of eqns. (fTT|) or (|12p by themselves are sufficient 
to model the propagation of the electric and magnetic 
fields. However, there are many traditional simplifica- 
tions which can be applied, and which in other treatments 
are even sometimes required in order obtain a simple evo- 
lution equation. In particular, the various envelope equa- 
tions all use co-moving and/or envelopes as a 
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preparation for discarding inconvenient derivatives: Here 
such steps are optional extras. 

These are all considered in more detail for a factorised 
wave equation approach in Ts'] , but here I have adapted 
them for this context. 



A co-moving frame can now be added, using t' ~ 
t— z/vf. This is a simple linear process that causes 
no extra complications; the leading right-hand side 
(RHS) iar(3rUJ = ikr term is replaced by i[arj3rUi =F 
kf), for frame speed Vf = cuo/kf. Setting kf — 
kr(j-tJo) will cancel the phase velocity Vp of the pulse 
at Wo, not the group velocity. 



D. Diffraction 

One important feature lacking in this approach is the 
handling of transverse effects such as diffraction, al- 
though they can be inserted by hand (at least in the 
paraxial limit) by adding the term %{d'^ + dy)F^/2kr to 
the RHS 12]. However, no treatment of transverse effects 
has been achieved in a native directional fields descrip- 
tion on the basis of the first-order equations - although 
transverse terms arise naturally in the second-order equa- 
tion resulting from taking the curl of eqn. Treating 
nonlinear diffraction [sl] suffers the same difficulties, al- 
though presumably it might be incorporated in an anal- 
ogous way as to ordinary diffraction. 



2. The field can be split up into pieces localized at 
certain frequencies, as done in descriptions of opti- 
cal parametric amplifiers or Raman combs (as in, 
e-g-, @, in, nil)- The wave equation can then be 
separated into one equation for each piece, coupled 
by the appropriate frequency-matched polarization 
terms (see, e.g., [26j). 

3. A carrier- envelope description of the field can 
easily be implemented with the usual prescrip- 
tion of [13, F+[t) = A{t) eyi^[i{ujit - kxz)] + 
A*{t) exp[— i(a;it — kiz)] defining an envelope A{t) 
with respect to carrier frequency uji and wave vec- 
tor fci; this also provides a built-in a co- moving 
frame Vf = wi/fci. Multiple envelopes centred at 
different carrier frequencies and wave vectors (wi, 
ki) can also be used [26, ■28i]. 

4. Bandwidth restrictions might be added (see below), 
either to ensure a smooth envelope or to simplify 
the wave equations; in addition they might be used 
to separate out or neglect frequency mixing terms 
or harmonic generation. As it stands, no band- 
width restrictions were applied when deriving eqns. 
(fTTj) or - there are only the limitations intro- 
duced by the dispersion and/or polarization models 
to consider. Typically we would expand the model 
parameters to the first few orders about some con- 
venient reference frequency wq- 



5. Mode averaging is where the transverse extent of a 
propagating beam is not explicitly modeled, but is 
subsumed into a description of a transverse mode 
profile; as such it is typically applied to situations 
involving optical fibres or other waveguides. Thus 
we could use mode averaging when calculating the 
effective dispersion or nonlinear parameters. See, 
for example, [2^ for a recent approach, which goes 
beyond a simple addition of a frequency depen- 
dence to the "effective area" of the mode, and gen- 
eralizes the effective area concept itself. 



IV. THIRD-ORDER NONLINEARITY 

Third-order nonlinearities are common in many mate- 
rials, for example in the silica used to make optical fibres 
(see, e.g., [gI]). There are many applications of significant 
scientific interest, for example, white light supercontin- 
nua [3ll - l33j . optical rogue waves [13); or filamentation 

Here we study propagation through such a material, 
with non-reference linear responses k^, (describing e.g. 
loss and dispersion), and instantaneous magnetic third- 
order nonlinearity WA along with the more common 
electric type Xe- For such a system, and for plane- 
polarized fields, the propagation equation is 



-ikr 



ikr eo 



9[F+'{t)]*F^ 



(13) 



This is a generalized nonlinear Schrodinger (NLS) equa- 
tion, and it retains both the full field (i.e. uses no enve- 
lope description) and the full nonlinearity (i.e. includes 
third- harmonic generation). The only assumptions made 
are that of transverse fields and weak dispersive and non- 
linear responses; these latter assumptions allow us to de- 
couple the forward and backward wave equations. This 
decoupling allows us, without any extra approximation, 
to reduce our description to one of forward-only pulse 
propagation. The specific example chosen here is for a 
cubic nonlinearity, but it is easily generalized to the non- 
instantaneous case or even other scalar nonlinearities. 

We can transform eqn. (fT5|) into one closer to the 
ordinary NLS equation by representing the field in terms 
of an envelope and carrier 



F+{t) ^ A(t)exp [i {ujot 



koz)] + A* (t) cxp [—1 {ujot — koz)] , 

(14) 



where we choose the carrier wave vector to be fco ~ 
kr{ujo) — u!oar{ujQ)l3r{ujo). After separating into a pair of 
complex-conjugate equations (one for A and one for A*), 
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and ignoring the off-resonant third-harmonic generation 
term, this gives us the expected NLS equation without 
diffraction. The chosen carrier effectively moves us into a 
frame that freezes the carrier oscillations, but this phase 
velocity (vp = uj/k) frame differs from one that is co- 
moving with the pulse envelope (i.e., one moving at the 
group velocity Vg — dui/dk). After we transform into a 
frame co-moving with the group velocity at wq, where 
Ag — uji^lvg"^ [uj^)) — v~^{luq)\, the wave equation for A{uj) 



d,A = +iK(uj)A 



2x\A(t)\'A{t) , (15) 



where K{lo) — kr[n^{Lj) + Kp(tj)]/2 + A^; and x — 
Xe — (/xoe^/eo/ir)XA'- th&t has been assumed to derive 
this standard envelope NLS equation is uni-directional 
propagation and negligible third-harmonic generation. 
The self-steepening term, often seen in (or added to) 
NLS equations arises from the frequency dependence of 
kr- This self-steepening has both electric and magnetic 
contributions, which can be adjusted independently, as 
has been pointed out by Wen et al. Q for the case of 
the SVEA limit. In Section IVIl I discuss how the impor- 
tance of each contribution varies with frequency for both 
a double-plasmon model (as in Q), and a wire- array and 
split-ring model more typical of practical metamaterials. 

It is worth comparing this eqn. (jl5p to D'Aguanno et 
aVs eqn. (5) [hereafter eqn. (DMB5)]. Although in 
many respects they appear to be the same, mine is far 
more general and can be applied (at least in principle) 
to an arbitrarily wide pulse bandwidth, whereas theirs is 
subject to the rather restrictive SVEA. For example, my 
eqn. (fT3)) results from only one "slow evolution" approxi- 
mation, as opposed to the numerous steps, substitutions, 
and approximations in Section 2 of [sl. I also retain the 
possibility of arbitrary dispersion K{uj), whereas theirs 
retains only the second-order part (i.e. as cx , which in 
the frequency domain would be cx oj^). Indeed, with the 
dispersion and nonlinear factors in my eqn. (|13p com- 
bined, that full-spectrum wave propagation equation is 
scarecly more complicated than eqn. (DMB5). Similar 
remarks also hold when comparing eqn. (jl5p to Wen 
et al.'s [3]: but although Wen et al.'s result is also re- 
stricted by the SVEA, it does at least allow for diffrac- 
tion. Both, however, along with Scalora et aL's form Q, 
cannot model the full non-envelope field, nor revert to 
an exact and explicitly bi-directional form, as in my eqn. 
® or (dni). 



V. SECOND-ORDER NONLINEARITY 

Treating a second-order nonlinearity is more compli- 
cated than the third-order case, since it typically cou- 
ples the two possible polarization states of the field to- 
gether. Such interactions occur in materials used for op- 
tical parametric amplification, and have long been used 
for a wide variety of apphcations (see, e.g., |28l. [ssl [39|). 



To model the cross-coupling between the orthogonally 
polarised fields, it is necessary to solve for both field po- 
larizations; and to allow for the birefringence we need 
two pairs of (non- reference) linear responses, i.e. k^j,, 
and K'^x^ ^iiy- 

As an example, I choose a magnetic nonlinearity that 
couples Hy and in the same way as the electric non- 
linearity couples Ex and Ey^ although other configura- 
tions are possible. This means that the /J^u x P term 
in eqn. ([T]), which represents the non-reference part of 
the electric response, needs to include those for the stan- 
dard second-order nonlinear terms (here oc E^Ey 
and Py oc E'^). Similarly, the ar/^oM term has ones 
for the complementary second-order nonlinear magnetic 
response. Note that second-order nonlinear magnetic ef- 
fects have been measured in split ring resonators by Klein 
et al. (13, 113 . 

Since it is convenient, I split the vector form of the 
wave equation up into its transverse x and y components. 
By noting that the definition of means that iJ+ = 
—Fj'ar/ I3r, the ID wave equations can be written as 



F+ 



.ikr ep 
' 2 er 



fxo er 

£0 \Mr 



j[F+{t)F+{t)] 
(16) 



1 



Mr 



(17) 



These wave equations for the field are strikingly similar to 
the usual SVEA equations used to propagate narrowband 
pulses; the main differences are the addition of terms 
for magnetic dispersion (k^^, i^^iy) and nonlinearity (x^), 
and the lack of a co-moving frame. 

We can transform eqns. ([T51) and pT|) into a form 
close to the usual equations for a parametric amplifier by 
representing the x and y polarized fields in terms of three 
envelope and carrier pairs: 

F+(t) = ^i(i)exp[z {uJit- kiz)] -I- ^^(t) exp [-z (wit - kiz)] 

+ A2{t) exp [l {i02t - k2z)] + ^2(0 exp [-1 {i02t - k2z)] 

(18) 

F+{t) Asit) exp [i {uj3t - ksz)] + A;{t) exp [~i {uj^t ~ ksz)] , 

(19) 

where uj^ — uJi + uj2- After separating into pairs of 
complex-conjugate equations (one each for all Ai and 
A*), and ignoring the off- resonant polarization terms, we 
transform into a frame co-moving with the group velocity, 
although here we select the group velocity of a preferred 
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frequency component, with = uj{Vg^ — v^^). The 
wave equations for the Ai{uj) are then 

Zki 

(20) 

d,A, = +iK2{lo)A2 + '^x-^[^A^{t)Al{t)] e-^^'^^ 

2K2 

(21) 

= +^K:,{uJ)A:, + 'A^+y [A^(t)A2{t)] e+''^'=^ 

2K3 

(22) 

Here ifi,2(w) = ki,2[K^x{'^) + '«Ma(w)]/2 + Ag and 
K^{uj) = k^[K^y{Lu) + k;^2^(ci;)]/2 + Ag; we choose /c^ for 
each equation differently (i.e., with kr S {^1,^2,^3}); 
also the phase mismatch term is Afc — fcs — fc2 — 
ki. The combined nonlinear coefficient is = Xe =t 
{^la/eQ){er/^irf''^X^i■ 



VI. DISCUSSION 

Examining the respective roles of the reference permit- 
tivity tr and permeability in eqns. (fT3)) and (fT6)) . (IT7| . 
we see that as far as dispersion and other linear effects 
are concerned, the two components simply add. In con- 
trast, their effect on nonlinear terms is more dramatic: 
with the ratio — trliJ^r scaling the nonlinear correc- 
tions to the magnetization into the electric field units of 
F*. This is because determines how much of a given 
directional field is electric field and how much mag- 
netic field; large values of correspond to cases where 
the magnetic field is most prominent. Indeed, Y is just 
the reciprocal of the electromagnetic impedance of our 
chosen reference medium, and only if Y is real- valued do 
propagating fields exist, since otherwise the fields become 
evanescent. 

Figures [1] and [2] show how Y varies with frequency for 
two different metamaterial types, with the dispersions 
encoded on and \ir and scaled by uP' to moderate the 
low-frequency singularity of the Drude response. The 
extreme limits of large Y occur when ^ \^r\^ that 
is, usually just at an edge of a non-propagating band, 
where iir is about to change sign. In such a region, it 
would be better to revert to the fields, or to rescale 
the propagation equations into units of magnetic field 
(e.g., with some = 6^/2/3^). 

In previous work 0411 , a Drude type response for both 
e and /i was assumed, where tril^r oc 1 — ui1^^/{uj'^ — 
i^e.fii^)', and this situation is shown on Fig. [TJ However, 
although the dielectric response in metamaterials (e.g., 
a wire grid array [HI, E^l) often has this behaviour, the 
magnetic response of split ring resonators (SRRs) dif- 
fers. SRR magnetization is best described by a pseudo- 
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FIG. 1: Normalised electromagnetic field ratio and e, curves 
for Drude responses in both permittivity e,. and permeability 
^r- The magnetic resonance at o;^ — u>^/y/2 is lower than 
the (dielectric) plasma frequency cj^, and 7e = 7^t = O.Ola;^. 
Large \Y\ corresponds to mostly magnetic fields, small 
yj corresponds to mostly electric G*. The frequency ranges 
of propagating negative phase velocity (NPV) and positive 
phase velocity (PPV) light are shown. 




1 2 

CO/COg 



FIG. 2: Normalised electromagnetic field ratio and e, fj, curves 
for Drude response permittivity (dot-dashed line), with 
a pseudo-Lorentz response for the permeability (dashed 
line). The magnetic resonance at uj^ = uji/\/2 is lower than 
the (dielectric) plasma frequency u^, and 7^ = 7^/5 = Q.Qlujt, 
Ff_i = 5. The high-frequency behaviour of the pseudo-Lorentz 
model "illegally" increases faster than that of the Drude 
model, whereas a properly causal form should match it. Apart 
from detail, and the high-frequency behaviour, replacing the 
pseudo-Lorentz model with the (causal) Lorentz one gives a 
figure of similar appearance. The frequency ranges of propa- 
gating NPV and PPV light are shown. 

Lorentz modeP [43] with /i^ fx l + F^a;^/(a;^ — w^ — 17^0;), 
although sometimes a true Lorentz response is used in- 
stead, cx 1 -I- F^w^/(w^ — — i7pw). Note the differ- 
ence between the numerators in these latter two expres- 
sions: a frequency-dependent uP' versus a constant ma- 
terial parameter o;^. The pseudo-Lorentz model has an 
incorrect high-frequency behaviour, and so it is incom- 



Also known as the "F-model" 
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FIG. 3: Normalised refractive index n (solid line) and group 
velocity Vg = dk/duj (dashed line) for the system defined in 
fig. [2I for only those frequencies where the field propagates. 
By comparing with fig. [2] we can see that the edges of the 
NPV band are dominated by the magnetic field, whereas the 
lower edge of the PPV band is dominated by the electric field. 
Note that in the NPV band, the sign of the group velocity is 
not tied to the sign of the phase velocity. 



patible with the Kramers-Kronig relations that enforce 
causality. However, at low and medium frequency it is a 
better match to the physical response of SRRs, and so I 
use it for Figs. [2] and [31 

There are frequency ranges over which the linear ma- 
terial responses vary dramatically, and in particular on 
Fig. [3] (which uses the same model as Fig. [2]) this is 
evident for both the refractive index n and group veloc- 
ity Vg. If we aim to operate in such regions, this leads 
to two potential complications. First, if we have chosen 
reference parameters that do not match the linear ma- 
terial responses exactly, then the correction terms will 
become large, meaning that our wavelength-scale "slow 
evolution" approximation may come under threat. Sec- 
ond, even if our reference parameters do match the linear 
material responses exactly, our wavelength-scale will have 
become frequency dependent, and so again our "slow evo- 
lution" approximation may be threatened. In either case 
the solution is simple - we just need to revert to the bi- 
directional wave equations [i.e., eqns. ^ or (ITUl) ]. How- 
ever, this does not necessarily mean that any backward 
evolving fields are generated (as explained in [11| , and fol- 



lowing a different approach in [15|), so that in principle 
one could optimize the propagation by reinstating it only 
over those frequency ranges where it becomes necessary. 



VII. CONCLUSIONS 

I have derived a uni-directional optical pulse propaga- 
tion equation for media with both electric and magnetic 
responses, based on the directional fields approach 
This involved a re-expression of Maxwell's equations, and 
required only a single approximation to reduce a one di- 
mensional bi-directional model, to a uni-directional first- 
order wave equation. The simplicity of this approach 
makes it very convenient in waveguides, optical fibres, or 
other collinear situations. The important approximation 
is that the pulse evolves only slowly on the scale of a 
wavelength; and indeed this is a valid assumption in a 
wide variety of cases - note in particular that nonlinear 
effects have to be unrealizably strong to violate it [Tij. 
The result has no intrinsic bandwidth restrictions, makes 
no demands on the pulse profile, and does not require a 
co-moving frame - unlike other common types of deriva- 
tion [i-ll. 

The resulting equations have the advantage that they 
are straightforward to write down, despite containing the 
complications of both electric and magnetic responses, 
and that a carrier-envelope representation or co-moving 
frames are easy to apply if desired, requiring no fur- 
ther approximation. In this, they match the clarity 
and flexibility of factorized second-order wave equations 
[isl li^ . lis! , but they can more easily incorporate the ef- 
fects of magnetic material responses - albiet at the cost 
of being restricted to one dimensional propagation. 
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transform into frequency space: 

V X U{t) = +dter ★ E{t) + 3{t), 

V X E(t) = -dtfir * H(t) - MoK(t) (Al) 

V X H(a;) -luj ar{uj)^E{uj) + 3{oj), 

V X E(a;) = +tuj /3^(w)2H(w) - HqK{uj). (A2) 

I now rotate the V x H equation by taking the cross 
product with u, 



u X (V X H) = -zo; (u x E) + u x J, 



(A3) 



scale each part by /3r and respectively, while insisting 
that these parameters do not depend on position. Thus, 

U X (V X Prii) = -iOJ PrCtl (u X E) + U X J, 

V X a^E = +ILJ UrPlH - q;,,^oK, (A4) 
and then take the sum and difference - 
V X a^E ± u X (V X /?rH) 

= +ILO UrPlVL T PrCtl (u X E) 

± u X /3r J — arfJ-oK. (A5) 



The vector fields are defined in eqn. but that 
neglects the longitudinal part of H. Thus, for complete- 
ness, we also need to define G° = u • /3rH (as in (lTl.l46l|). 
Their form means that I need to convert both the sec- 
ond term on the LHS of the sum-and-difFcrence equation 
above, as well as the RHS. It is most important for the 
LHS to be simple, because this will define the type of 
propagation specified by the RHS. In the following I use 
the vector identity. 



u X (V X H) - V (u • H) = V X (u X H) , (A6) 

along with u x [u x H] = [u • H]u — H, so that 

u X = u X ar'E ± u X [u X /J^H] (A7) 

= u X ar'E =F /3rH ± [u • /3rB.] u (A8) 

= uxa^E=F/3rH ±uG°. (A9) 



Continuing the derivation. 



V X a^E ± u X (V X /3rH) = +iujarPl'H. =F lUJ/Sral (u x E) ± u x /J^ J - ar/LtoK, 
V X a^E ± V X (u X PrU) ± V (u • /3rH) = -l-iwar/S^H =f lOjfUral (u x E) ± u x /SrJ - q;,.^oK 



(AlO) 
(All) 

V X [a^E ± (u X /3,-H)] — +iLdar /3^H =F lOjPral (u X E) =F V (u • /33) ± 

(A12) 

V X G=F = {arPr^ T /3ra^ (u X E) } T V (u • ;3^H) ± u x J - a^/zoK 

(A13) 



I now rearrange eqn. (|A13|) to give the final form, in which I substitute dtP = J(t) and 9tM = K(<) to match eqn. 
(P). Thus 



X G"^ = luj {arl3ruG° - (u X [u X H]) =F /Sr^r (" x E)} T VG" ± u x /^^ J - a^^oK 



(A14) 



=F«w {/3rQ;^ (u X E) ± Qfr/Sr (u X [u X H])} +tujarl3ruG° =F VG° ± u x /^^ J - a^MoK, (A15) 



and finally 



V X ±iLo a,./?,.u X G^ + tujar(3ruG° ± VG° =f u x /SrJ — ar/xoK 

= ±tUJ arPrVi X G^ + lUJarf3rUG° ± VG° ± lOj/SrU X P -I- lWQ!r/ioM. 



(A16) 
(A17) 



Note that for J{t) = dtP = dtKgE{t), where is some 
complicated but scalar dielectric response function, we 



have 

=FU X (3r3{uj) — ±lUjl3rU X P = ±lU}al(3rHe ★ U X E 

(A18) 



= ±^uj^K,*{ux [G+ + G-]) 



(A19) 

and for K(i) = i9tM — 9tK^H(t), where is some com- 
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plicated but scalar magnetic response function, we have 

— arA*oK(aj) = +iiOfioar'M. = +2UjfiQarKf^ -kH (A20) 

= +iuJHQarK^ * (u [u ■ H] — u X u X H) 

(A21) 

ar 

(A22) 

^tLu^io^K^ * (u X [G+ - G-] - 2uG° 

(A23) 



+iaJ/J,0-r^Kp * (uG° — U X [u X /3.rH] 



Finally, when generating eqn. (jASp . we lost the longi- 
tudinal part of V X H = 9erE + J (i.e. that parallel to 
u). This is 

u • V X H = -lujalu • E + u • J 

(A24) 

u • V X (G+ - Q- + 2uG°) = -lujarPrU ■ (G+ - G") 

+ 2/3rU • J (A25) 
2u • (VG° - u X VG°) = -iLuar/SrU ■ (G+ - G") 

+ 2/3^u • J (A26) 

2U • VG° = -lOJUrPrU ■ (G+ - G") + 2PrU ■ J, (A27) 
since V x G°u = G° V x u - u x VG° = -u x VG°, and 

u • V X (G+ - G") = lujarPrU ■ u X (G+ + G") 

+ 2u • VG° + 2iujf3rU • u X P 
(A28) 

= 2u-VG°. (A29) 

Appendix B: Correction terms 

In this appendix I work through the details of how 
the polarization and magnetization terms scale with re- 
spect to one another. To simplify matters, I assume all 
corrections are scalar since when and fir are not field- 
polarization or orientation sensitive, the scalings remain 
the same, even if the specific field terms may vary (e.g. 
ExEy instead of E"^). 



Consider the general unidirectional equation for 
(i.e. eqn. (fTU)) ). and replace the polarization and mag- 
netization terms with dimensionless response parameters 
Qe and multiplied by the appropriate field E^- or Hy. 
Then replace E^ and Hy with their representation in 
terms of , so that 



2 ar 



ILO ILO Pr 

Px + —fiQMy = — — q^eQE^ + —qf_,fiaHy 



lUJ 



(Bl) 



y 



1 17' + 

Ur Pr 



lUJUrPr 



ikr Co 

~2~Z 



9e£0 + 



(B2) 

P^r " 

(B3) 



9e +9^ 



Mo £7- 
eo fir 



F+ 



(B4) 



remembering that F''" = E = (/3r/ar)H, and that — 

ttr, and fir = Pr- 

Since we consider the electric- field-like field F+, the 
polarization corrections are trivial to write down; as for 
an m-th order nonlinear term, g,; = Xji^"'"*^™"^^ This 
means we need only concentrate on the magnetization 
correction. If is that for an m-th order nonlinear term, 
then qf, = Xt^H''-^ = xMr / Prf"^'^^ F+^'^-^'K Writing 
down only the term in square brackets from eqn. (jB4p 



gives us 



eo fir \ Pr 



F^ 



(B5) 



A*0 / 



£0 V ft-r 



+ (m-l) 



F^ 



(B6) 



Note that corrections for linear loss or gain are first-order 
processes (i.e. with m = 1), where for loss we need q ^ 
«7, with 7 > 0; Thus for loss the whole correction term 
will be proportional to — 7F+, as would be expected. 
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